Graphene in periodic deformation fields: dielectric screening and plasmons 
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We consider the effect of periodic scalar and vector potentials generated by periodic deformations 
of the graphene crystal lattice, on the energy spectrum of electrons. The dependence of electron 
velocity near the Dirac point on the periodic perturbations of different types is discussed. We also 
investigated the effect of screening of the scalar potential by calculating the dielectric function as 
a function of the wave length of the periodic potential. This calculation shows that the periodic 
scalar field is strongly suppressed by the screening. Using the dependence of electron velocity on 
the periodic field we also studied the variation of the plasmon spectra in graphene. We found that 
the spectrum of plasmon excitations can be effectively controlled by the periodic strain field. 
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I. INTRODUCTION 

The enormous interest to graphene is related to 
the unique physical properties of this two-dimensional 
materialir^ which most possibly will be used in the fu- 
ture in numerous technological applications^—. An ex- 
ample of such properties of interest for the applications 
is the existence of very unusual spectrum of plasmon ex- 
citations with THz frequencies^, that can be used in 
optoelectronics and communications. 

One of the most important problems to be solved to 
use graphene in electronics is the realization of effective 
control of the parameters of energy spectrum such as the 
electron energy gap and/or the velocity of electron and 
holes. It is already known that by using the electrostatic 
gating one can vary the carrier density of graphene (i.e., 
the location of chemical potential) . Recently, it was also 
proposed to use the external strain field to change the 
energy spectrum - this is called the strain engineering of 
graphene^iISr— . The idea is mostly based on unusually 
strong effect of the external deformation acting on the 
energy spectrum of graphene quite similar to the external 
electric and magnetic fields 3 -^. 

In this work we consider the effect of periodic fields, 
which can be generated by periodic deformations, on the 
energy spectrum, screening of electron-electron interac- 
tion, and on the plasmon excitations in graphene. It 
was already pointed out that such periodic fields do not 
open the gap near the Dirac point but affect the veloc- 
ity parameter of electrons and holes^— making the en- 
ergy spectrum anisotropic. Recently, the effect of peri- 
odic modulation on electron spectrum of graphene has 
been studied experimentally-^. Here we reconsider this 



problem in more details, concentrating on possible coex- 
istence of periodic scalar and vector potentials, which, to 
our knowledge, has not been done before. This is im- 
portant since a generic deformation produces both scalar 
and vector potentials 3 - : 12 ' 13 . We use a different method 
to solve the problem^, and find that our numerical re- 
sults are in agreement with those of Refs. |14Hl6| . We 
also discuss the role of screening due to the free carriers 
in graphene, and we find that screening can substantially 
suppress the effect of periodic scalar potential. 

The combined effect of periodic scalar potential and 
constant magnetic field in graphene has been recently 
investigated by Wu et al J£ They found that the structure 
of Landau levels can be also strongly affected by the one- 
dimensional (ID) periodic fields. 

The physics of plasmons in graphene has been inten- 
sively discussed recently by many authorsi 2 ^— The cal- 
culations have been performed in frame of standard RPA 
approximation taking into account the electron energy 
structure of graphene near the Dirac points, as well as for 
the whole energy spectrum at the honeycomb lattice 2 ^. 
The effects of magnetic field, finite temperature, chemical 
potential have been investigated in the same approach 26 . 
In this work we discuss the effect of external periodic 
fields on the screening and on the energy spectrum of 
plasmons. 



II. ELECTRON ENERGY SPECTRUM IN 
PERIODIC FIELDS 

We consider first the transformation of electron en- 
ergy spectrum related to deformations of the graphene 
lattice. It is known that the deformation of graphene 
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is equivalent to the generation of electric and magnetic 
fields, which can be described by scalar V(r) and vector 
A(r) potentials 3 ' 13 . The relations between the compo- 
nents of strain tensor t%(r) and the scalar and vector 
gauge potentials ar o 27 ' 28 



V(r) = g (u xx + u yy ), 

A x (r) = — (u x: 
a Q 



2f3t 

u yy)i A y (r) = u xy , (1) 

a 

where g is the deformation potential, t is the hopping en- 
ergy, the parameter (3 is defined by (3 — — <91ogi/<91oga , 
and oo is the lattice constant. In the following we con- 
sider one-dimensional periodicity of the deformation, and 
in view of Eq.(l) we assume that the periodic in x strain 
field generates periodic scalar and vector potentials V(x) 
and A(x). 

The Hamiltonian of electrons in graphene near the K, 
Dirac point in periodic scalar and vector fields reads 



H = -ivcr ■ (V - iA) + V, 



(2) 



where the Pauli matrices er act on the sublattice label 
and we use the units % = e = 1. This Hamiltonian 
describes low-energy excitations of the electronic system 
in graphene. 

The corresponding Schrodinger equation for spinor 
wave function ip ( r ) — (<p, x) is 



V 

- vA 



+ 



ivd- 
s - 



- vA- 
V 



= 0, 



(3) 



where d± = d x ± id y and A± = A x ± iA y . Since the 
potentials A(a;) and V(x) do not depend on y and depend 
periodically on x, we take ip,\ ~ e . 

In frame of the k ■ p approximation 2 ^, we have to cal- 
culate first the wavefunction i/j(x) at k = 0. The corre- 
sponding equations for the spinor components at k = 
are 



(e — V)ip + ivx + vA-x z 
ivtp' + vA + tp + (e — V)x - 

From Eq. (5) follows 

iv(p' vA + ip 



0, 
= 0. 



X = 



e-V 



V 



(A) 

(5) 



(6) 



where prime means d x . Substituting Eq.(6) into Eq.(4) 
we obtain the following equation for (p(x) 



(e - V)ip + 



V 2 ip" 



v 2 V'tp' iv 2 A' + ip iv 2 A + tp' 



e-V (e-V) 2 e-V e-V 
iv 2 V'A + tp iv 2 A-tp' v 2 A + A-(p 



(e-V) 2 e-V e-V °' (7) 

Let us assume the existence of solutions of Eq. (7) with 
= (here we assume that V ^ 0). Then we get 

iV'A-L. 



" ( y + ->i.\ 



V 



-A + A_)tp = 0. (8) 



This is the equation for ip(x) in the AC Dirac point, cor- 
responding to the lowest energy band. 

Let us consider first some particular cases, when the 
periodic field is purely scalar V(x) or purely vector field 
A(x). 



A. Periodic scalar potential 

In the case when A(x) = and V(x) ^ 0, Eq.(8) 
essentially simplifies to 



and has two different solutions 



<Pi,2(x) = exp ( ±— / V(x')dx 



Correspondingly, using Eqs.(6) and (10) we obtain 



XiA x ) = Tcxp ±- / V(x')dx' 



(9) 



(10) 



(11) 



Then the normalized basis functions in k ■ p approxima- 
tion are 

Mr) = ^(^), Mr) = ^(Zl), (12) 



2S \ Xi 



25 V X2 



where S = C x C y is the area of graphene sample. Calcu- 
lating the matrix elements of the Hamiltonian (2) with 
basis functions (12) we find the effective Hamiltonian 

\ v(^ik x - j 2 k y ) vk x , )' 

where we denote 

71 = 1 J dx cos J V(x') dx'^j , (14) 

72 = \ J dx Sin (~ / V ^ dx '^j ' ^ 

and L is the period of the potential V(x). In correspon- 
dance with Eqs.(14) and (15) both parameters 71, 72 < 1. 

The Hamiltonian (13) describes low-energy spectrum 
in the periodic scalar field. It has the following eigenval- 
ues 



£1,2 (k) = ±v J k 2 + (jikx - i2k y Y 



(16) 



Taking on alternate k y = and k x = we find that due 
to the periodic scalar field, the components of electron 
velocity in directio ns x an d y, are renormalized, respec- 
tively, as v x /v = VI + 7J and v y /v — |7 2 |. Thus, in this 
case we always obtain v x > v and v y < v. 
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FIG. 1: Dependence of the renormalized velocity factor v/v 
on the amplitude Ao of periodic potential A y (x) at different 
values of periodicity parameter L. 



B. Periodic vector potential 

In the case when V(x) = and A(x) ^ we can use 
Eqs.(4) and (5) to find directly from these equations the 
spinor components of the wave function in the K, point, 
k = 



<p(x) = exp I i j A + (x ) dx j , 

X(x) = exp (i A- (V) dx'^j . (17) 

Then we can introduce two k • p basis functions in the 
form 



^(r) = iV ie ik - r (J ) , Mr) = A 2 e 4k ' r ( 



(18) 



where Nj arc the normalization factors 



C y / exp ( =f2 / A y (x') dx' I dx 



1 -1/2 



(19) 



for the crystal of size C x x C y . Calculating the matrix 
elements of the Hamiltonian (2) with V(x) = in the 
basis of k -p functions (18) we get 



H 



vk_ 
v*k + 



where 



v = vNiN^Cy I ip* (x) x(x) dx . 
Jo 



(20) 



(21) 



Using Eqs.(17) and (21) we finally obtain 



L 



X 



exp 



2 J A y (x')dx'\ dx x 

X 2 \ 1-1/2 



Ay(x') dx' I dx2 



o 



(22) 



The dependence of v/v from the amplitude A Q of the peri- 
odic potential A y (x) — A$ sin(27rx/L) at different values 
of period L is presented in Fig. 1. The renormalized elec- 
tron velocity decreases in the periodic vector field. This 
is in agreement^ with Ref. [Ty] . 



C. General case: both scalar and vector potentials 
are nonzero 

In the general case when both V(x) ^0 and A (a;) 7^ 
we cannot find simple analytic solutions but we can 
analyze further Eq.(8) for f(x). For this purpose we 
present this equation as 



<p" + a<p' + bip = 0, 



where we denoted 



a(x) 



YL 
V 



2iA, 



V 2 

b(x) = - iA' + 
v^ 

Then after substitution 



iV'A+ 
V 



-A+A- 



ip(x) = f(x) exp f -- J a(x') dx' 
we obtain the equation for the function f[x) 



(23) 

(24) 
(25) 

(26) 



/ = o. 



(27) 



2 J V 4 8 2/ 
This is the Schrodinger equation for a particle of unit 
mass with energy e = in the potential 



V" 



3(V') 2 
8V 2 



A y 



A' V'A„ 



2V 



V 2 

*?• (28) 



Here we note that U(x) does not depend on A x . It can 
be used when A y = since in this case we can take the 
solution for fix) corresponding to ip{x) from Eq. (10). 
This way we obtain simple generalization of Eqs. (10) 
and (11) for A x ^ and A y = 



fi,2(x) = exp 
Xi,2( x ) = Texp 



± 



V(x') 



± 



+ A x (x')j dx' 
l + A x (x')J dx' 



(29) 
(30) 



Turning back to Schrodinger equation (27) we also note 
that the potential (28) is real and periodic, U(x) = U(x+ 
L), so that we can present it as 



N 

U{x) = E ( 



(31) 
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where the coefficients u n can be found from the specific 
shape of potential U(x). 

We are looking for a periodic solution f(x) of Eq.(27), 
which can be presented as 



27rimx / L 



(32) 



with to is integer. Then using Eqs.(3f), (32) and (27) we 
find the matrix equation for the coefficients f m in Eq.(32) 



^ ^ Anm fm 0; 



(33) 



where A nm = \ k*5 nm + u„_ m + u* m _ n with k n = 2irn/L 
and u n = for any n < 1. 

One can also look for the solution of Eq.(8) in the form 
tp(x) = e ls ^ x \ This representation can be more conve- 
nient for numerical calculations with arbitrary periodic 
functions V(x) and A(x). In this approach we get the 
first order differential equation for £(x) 



V 



2iA x £ 



iV 2 



A' 



V'A^ 
V 



-iA+A- = 0, (34) 



where £(x) = s'(x). Note that the transition to the case 
V(x) —> formally corresponds to V'/V —> oo in Eq.(34). 

One can assume that like Eqs.(12) and (18), in the 
general case there are also two solutions of Eq.(8) for the 
envelope function ip(x). Then there are also two different 
solutions of Eq.(34), £i(x) and Correspondingly 
we get two solutions for the first spinor component 



<Pi(x) = cxp \i J £i(x)dxj, i = l,2. 

Then using Eq.(6) with e — we find the other com- 
ponents Xi( x ) an d X2{x). The obtained spinor function 
(tfiii Xi) T should be properly normalized. As before, we 
use these independent solutions for our k ■ p basis pre- 
sented by Vi(r) = Nie lU ' T (vi,Xi) T - 

Thus, in the general case of arbitrary periodic pertur- 
bations we obtain the effective Hamiltonian 



, 2va!k x + 2va 2 k y v-yk- + v5k + , . 
1 ~ 1 vj*k + + v5*k- 2vfi 1 k x + 2v(3 2 k y ) ' 1 > 



where we denote 



a\ + ia 2 = N 2 C y / (p* xi dx 



p 1 + ip2 = NlC y <p* 2X 2dx, 



C = N^Cy / <p\ X 2 dx, 



(36) 



5 = NiN 2 C y / xl </?2 dx. 



> 

>> 
■> 
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FIG. 2: Dependence of the renormalized velocity factor v y /v 
on the amplitude of periodic scalar potential V(x) at different 
values of periodicity parameter L. 



The eigenvalues of Hamiltonian (35) are 

£1,2 (k) = v(ai + 0i)k x + v(a 2 + (3 2 )k y 

±v[( ai - P.fkl + (a 2 - p 2 fk 2 y + (|C| 2 + \S\ 2 )k 2 

+2Re(C<T)(fc 2 ^ k 2 y ) +4Im((6*)k x k y ] 1/2 . (37) 
For k x = we obtain £i, 2 (k y ) = v\ 2 k y , where 

tifl/v = («2 + fa) ± [(a 2 - (3 2 f + (|C| 2 + \S\ 2 ) 

-2 Re (05*)] 1/2 , (38) 

and for k y = we get £\, 2 {k x ) = vf 2 k x , where 

vfjv = (a, + ± [(at - (3 X ) 2 + (\(\ 2 + \S\ 2 ) 

+2Re(C5*)] V2 . (39) 

Using Eqs.(36), (38), and (39) one finds the components 
of electron velocity in graphene in the case of arbitrary 
periodic perturbation described by the fields V(x) and 
A(x). 

D. Longitudinal standing wave 

For a longitudinal strain wave, the components of de- 
formation tire Ux — 

(x) and u y — 0. Then in accor- 
dance with Eq.(l) we get A v = and due to Eq.(28) we 
can use the solutions (29) and (30). 
Using Eq.(36) we find 

a =-1/2, = 1/2, 

C = ~ S = 2Z I ^t'vJo V{ - X>) ^ (40) 
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and it follows from Eqs.(38) and (39) that v x = v and 
Vy/v = 2|C|. The dependence of v y /v on the amplitude 
of the periodic potential V(x) = Vq s'm(2irx / L) is shown 
in Fig. 2 for different values of periodicity parameter L. 



E. Transversal standing wave 

In case of transversal wave the only nonzero component 
of deformation is u y (x). Then in accordance with Eq.(l) 
we get V(x) — and A x {x) — 0. Correspondingly, the 
solutions for the components at k = are 

pi(ss) = exp(-^ A y (x')dx'y xi =0, (41) 

f2=0, X 2{x) = ex p(^ A y (x')dx'y (42) 

The solution for this case was already presented in 
Sec. 2B, see Eq.(22). 

It should be noted that in reality the longitudinal and 
transverse phonon modes in graphene are not completely 
independent - there is some mixing between them2£. 
The above consideration is fully justified for the phonon 
(standing waves) with small q. 



III. SCREENING 

The scalar potential V{x) generated by deformation 
wave in graphene is screened by electrons and holes. 
Screening is the main many-particle correction to the 
bare potential which should be taken into account. 

For this purpose, using the RPA approximation, we 
calculate the loop diagram presenting the polarization 
operator— 



. . m f cPkcfe „ 

n (q) = -iTry -^s-Glk 



where the Green function 



q, £ )G(k, £ ), (43) 



G(k,e) 



k 



(e + fj. + id sgn e) 2 — e\ 



(44) 



corresponds to the Hamiltonian of graphene without any 
perturbations, fi is the chemical potential, and £fc = vk. 
For definiteness we assume \i > 0. This quantity has 
been calculated in many papers (see, e.g., Refs. [20I32I — 
Hil) but we present here some intermediate expressions 
to discuss a generalization to the anisotropic case. 

Substituting Eq. (44) into Eq. (43) and integrating over 
e we find 



n (q) 



(27T) 



£k+q) - /(Ek)] 



-k+q 



v 2 k • q 



£k+q( £ k+q - £ V) 



-[1 - /(Ek)] 



-E k+q £ k + E 2 . + v 2 k ■ q 



£k£k+q(£k+q + Ek) 



(45) 
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FIG. 3: Dependence of dielectric constant e on q for different 
values of chemical potential fi. 



where f(e) = { exp[(e — fi)/T] + 1} is the Fermi dis- 
tribution function. In the limit of small 5 « (i/u we 
get 



n (q) 



d 2 k /(g k+q ) - /(e k ) 

(27r) 2 £ k +q - Ek 



(46) 



which gives us Ho(q — > 0) = — where v(e) = e/2irv 2 
is the density of electron states with energy e counted 
from the Dirac point. The second integral in Eq.(45) 
includes vacuum screening because it is nonzero at jj, = 0, 
i.e., in graphene without any carriers. In the limit of 
q — > this contribution to Ilo(q) disappears. 

In frame of the RPA, the screened potential V(q) is re- 
lated to the bare potential Vo(q) by V(q) — Vo + VoHqu = 
Vq/s, where uo(q) and u(q) refer respectively to bare 
and renormalized Coulomb interaction, and e(q) = 1 — 
uo(q)Ho{q). The dielectric constant e(q) calculated nu- 
merically using Eq.(45) is presented in Fig. 3. As we see 
in this figure, the scalar potential V(q) is substantially 
suppressed by the screening since e(q) 3> 1. 

The same dielectric function e(q) determines (within 
the RPA) screening of the electron-electron interaction, 
u(q) = uq/b. The value of the effective screening radius 
is R c = -(27re 2 n (0)) _1 = (27re 2 ^(^)) _1 . 

In the periodic field, the components of velocity are 
renormalized, so that the density of states changes from 
v(e) — e/(2TTv 2 ) to v(e) — e/(2-Kv x v y ). Taking into ac- 
count a variation of the chemical potential in the periodic 
field at a constant density of free carriers, ft, = \x^Jv x v y jv 
we find for the density of states at the Fermi level 
v = [i I {2'Kv^Jv x v y ). It means that v grows with decreas- 
ing carrier velocity as v ~ vj \Jv x v y leading to decreas- 
ing screening radius R c ~ ^Jv x v y /v. In other words, the 
periodic potential effectively enhances screening of the 
Coulomb interaction between electrons in graphene. 
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The periodic-field-induced variation of e-e interaction 
can affect the many-particle renormalization of the Fermi 
velocity^— due to modification of the screening. As 
follows from presented above estimations, this effect is 
especially important when v/ sjv x v y 3> 1. 

On the other hand, the anisotropy of velocity generated 
by the periodic field can essentially modify the renormal- 
ization group (RG) equations of Refs. [3al37j . Indeed, by 
calculating the Fock self-energy diagram in the case of 
anisotropic spectrum with v x ^ v y we find 



4-7T 



|k-q|>fc F 



(fq v x a x {k x - q x ) + Vy<Jy(ky 

q 



where kp is the Fermi wave vector, ik = (vlk^+Vyk 2 ) 1 ^ 2 . 
Correspondingly, the velocity correction is Svi = 5v^' + 
5v, , where 



5v 



r,(i) 



Sv 



,-,( 2 ) 



e 2 v t 
An 



e 2 Vj 
Auk 2 



d 2 q 1 



k-q|>fc F 1 £ k- 



k— q|>fc F 

In Eq. (47) we take the limit k -» 

< 2 r, /' ,/,/ ' 1: 
4tt J kp q 



d 2 q k ■ q 

q 



5v 



r,(i) 



dd 



v 2 sin 2 i 



(47) 



(48) 



(49) 



whereas in (48) we have to take first k = (fc,0) for Sv x 
and k = (0, k) for 5v v , respectively, and after that take 
the limit of k —¥ 



e 2 v x 
An 



dq 



■2 k 



v 2 x cos 2 Ode 



fr q Jo {v 2 cos 2 9 + v 2 sin 2 6>) 3/2 : 

2tt 52™ 2 (1^ 



9 



w 2 sin - 



o (v 2 cos 2 + vl sin 2 0) 3/2 



(50) 



The integrals over q in (49), (50) run from kp to g ma x — 
1/ao- Since the field-induced renormalization of velocity 
refers only to region of small q < 1/L, we will divide each 
of these integrals to the part from kp to 1/L (where, as 
we found, the spectrum is anisotropic), and to the part 
from 1/L to 1/ao with v x = v y = v. 

Let us assume for definiteness that for the bare values 
ix/iy < 1. Using (49), (50) and dividing each of integrals 
over q in two parts we find the following many-particle 
corrections to the velocity 



Sv x 



e 2 £o , e 2 v x £ 



7T (' 



SI 



K(m) - R(ro) 



v 



^-T^ + v^H-PHl, (5i) 

where we denoted £o = l°g(i/ao), £ = 
log(l/fc F L), K(m) = J* /2 (l ~ m sin 2 e)- 1 / 2 ^ 
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FIG. 4: The characteristics of RG Eqs. (52). In the limit of 
v x — > 0, all of them go to oo. The renormalization due to 
the e-e interaction shifts the initial values of v x + e 2 £o/4 and 
v y + e 2 £o/4 along a certain characteristic to the right. 



(complete elliptic 



P(m) 
Jo 



integral 

,2 n 



(I — to sin i 



to sin 



the 

,2 , 



1st kind 38 ), 
>d6, R(m) = 



l-«x/««- As 



of 

)- 3 / 2 sin 

)~ 3 / 2 cos 2 #<i6>, and to 
we are interested in the limit of small kp, the integrals 
over q are calculated with the logarithmic precision 
asssuming kp <C 1/L<C 1/ao- 

The first terms in the right hand sides of (51) lead to a 



constant shift of the bare values & 



+ e 2 £ /4. Note 



that in the limit of v x = v y , Eqs. (51) coincide with the 
ones from Refs. (35ll37| . 

Using (51) we obtain the following RG equations 



dv x 

dVy 



e 2 v x £; 



K(to) - ^| R(to) 



e 2 i 



[K(m) - P(m)] 



(52) 



The corresponding characteristics in the sector < v x < 
v y of (v x ,v y ) plane are presented in Fig. 4. As we see, 
the e-e-interaction-induced renormalization leads to the 
effective isotropization of the energy spectrum, which has 
been broken by the periodic field. 



IV. PLASMONS 

First we calculate the real part of polarization operator 
IIo(q, uS) for to ^ assuming that Sk — vk. It allows 
to consider the plasmons in graphene (corresponding to 
the poles of dielectric function e(w,q)) without external 
perturbations, 

Ren (q,^) = -y ^{[l-/( £k+q )] 
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FIG. 5: Plasmon spectrum ui P (q) of graphene for different 
values of the chemical potential /x. 



£k+q(£k+q 



v 2 \t • q 



4] 



£k+q[(£k+q - W) 2 - c kJ 

£kW + v 2 k ■ q 



-[l-/(ek)] 24 



£k[(£k 



IT' 



(53) 



The dielectric function e(q, w) = 1 — = 22 - ReIIo(q, w) 
can be found using polarization operator (53). The plas- 
mon spectrum is calculated by solving numerically equa- 
tion e(q,cdp) = 0. It is presented in Fig. 5. At small 
q -C /^A>the spectrum is w p (g) ~ y/q in agreement with 
Refs. |2l[l39|. As we see from Fig. 5, at larger q the 
plasmon dispersion is linear with q. When fi — > 0, the 
plasmon spectrum is linear. 

As we demonstrated in Sec. II, the energy spectrum 
of low-energy excitations in graphene under the periodic 
perturbation can be described by an effective Hamilto- 
nian (35). After some unitary transformation T this 
Hamiltonian can be reduced to the form similar of non- 
perturbed grapnene 



T X HT = v x a x k x + v y a y k y . 



(54) 



with renormalized values of electron velocity. For exam- 
ple, in the case of longitudinal standing wave, we have 
T = e~ t7TrTy / 4 , v x = v and v y — 2v\Q with ( denned by 
Eq. (40). 

The polarization operator ReII(q, uj) in the periodic 
field can be found using the same formula (53) after scal- 
ing transformation fej = (v/vi)ki and qi — (v/vi)qi. Then 
we find 



ReII(q, u) — Reil (q, w). 

V X Vy 



(55) 



Since the renormalization of electron velocities is differ- 
ent for longitudinal and transversal waves, the plasmon 
spectrum is different in these cases, too. 



FIG. 6: Plasmon spectrum ui P {q) of graphene under periodic 
perturbation for different values of the electron velocity (with- 
out anisotropy like in the case of periodic vector potential). 
Here the chemical potential = 10 meV. 



As shown before, the plasmon spectrum at q — > is 
proportional to ^fq. It corresponds to ReIIo(g,w) ~ 
q 2 juS 2 . Then the scaling transformation presented above 
does not change the polarization operator at q — > 0. 
Thus, the variation of electron velocity does not affect 
the plasmon spectrum at small q. 

In Fig. 6 we present the results of numerical calculation 
of the plasmon spectrum for different values of renormal- 
ized velocity v/v. It corresponds, e.g., to the presence of 
periodic field A y (x). This figure demonstrates that only 
the linear part of the spectrum can be strongly affected 
by the periodic field. 

In the case of longitudinal wave and in the limit of 
q — > we obtain 



ReII(g,w)~CM | : ? + — 



ml 



(56) 



It leads to the anisotropy of plasmonic spectrum. 



V. ROLE OF INTERVALLEY TRANSITIONS 

It should be pointed out that our consideration of the 
plasmonic spectrum cannot be extended to q of the order 
of the vector of inverse lattice. The point is that when one 
starts from the tight-binding approximation to describe 
the electronic structure of graphene, the Hamiltonian of 
the Coulomb interaction has the following form (index a 
labels sublattices A and B) 

Hint = ^ c R.<r CR - u ° _ ) C R' , Cr 4, 

= c L c k-q, CT wo(q)4, CT ,c k ' +qi(TS (57) 

kk'qcrcr' 
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where k and k' are any points in the Brillouin zone. Con- 
sidering the states near the Dirac points JC and JC' we can 
present (57) in a different form 

H int = [ c L 4 c k-q,^wo(q)4, ., i c k - +q , ff0 - 

+ c L; c k-q,<Tj «o(q - Q) ^a'jCk'+w'i] , (58) 

where i,j labels the valleys, k and k' are measured from 
the corresponding Dirac points, and Q is the vector be- 
tween the points JC and JC' . This expression shows that 
for plasmon excitations with momentum of the order of 
Q, intervalley transitions should be taken into account, 
as was pointed out in Ref. [22j ]. 

VI. CONCLUSION 

We considered the variation of electron energy spec- 
trum near the Dirac point in graphene under periodic 
perturbation related to the scalar and vector gauge fields, 
which can be generated by the periodic deformations. 
The possible source of such deformation fields is a peri- 
odic strain wave like in the case of the ultrasonic wave 
in solid. We found that in the general case with both 
V and A / there exist solution for the renormal- 



ized electron velocity, corresponding to the anisotropy of 
the spectrum. The problem substantially simplifies in 
some particular cases. Namely, for pure longitudinal and 
pure transverse periodic excitations the solutions have 
simple form. 

We also considered the screening of the scalar potential 
and found that it is strongly suppressed, especially at 
small q. It means that the main perturbation affecting 
the electron velocity is the vector potential A. 

We calculated the plasmon spectrum of collective ex- 
citations in graphene in the presence of periodic exci- 
tations. We found that the plasmon spectrum can be 
strongly affected by the periodic field. For the longitu- 
dinal excitations, one appears the anisotropy of plasmon 
spectrum. 
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